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It is argued that the title of this paper represents a mis- 
conception. Contrary to widespread beliefs it is electromag- 
netic field modes that are "systems" and can be entangled, 
not photons. The amount of entanglement in a given state 
is shown to depend on redefinitions of the modes; we calcu- 
late the minimum and maximum over all such redefinitions 
for several examples. 



I. INTRODUCTION 

If we have two quantum systems, then we all know 
that their joint quantum state may be entangled. This 
is the case when there is no way to write that state as a 
product state (or a mixture of product states) of the two 
systems. The systems may be particles such as electrons, 
atoms, or ions. In the quantum theory of light, which 
is a second-quantized theory it is not the "particles", 
the photons, that form systems, but rather the EM field 
modes. A (pure) quantum state of a mode is an arbitrary 
superposition 



between modes but not between photons. If we have a 
state 



I*) 



n=0 



(1) 



with the integer n the number of photons in that mode 
and a n arbitrary complex numbers only restricted by the 
normalization condition. 

Two field modes A and B are described by spatially 
dependent (vector) mode functions Fj (r) and by creation 
and annihilation operators at, a,j for j = A, B jj]. These 
field modes are independent systems if the mode func- 
tions are orthogonal, 



dfF*(r)-F j (r) = S ij . 



(2) 



The mode operators then satisfy the standard bosonic 
commutation relations 



[a,i, at] = Sij, 

and the number states \n)j of mode j are given by 

(at)" 



(3) 



(4) 



Before moving on to more interesting states, we first con- 
sider a simple example that shows that entanglement is 



|V) = (|o)a|i)b + |i)a|o) b )/V2, 



(5) 



there is one ebit of entanglement |!| between the modes 
A and B, although there is only one photon. In particu- 
lar, when modes A and B are located in different spatial 
regions then the entanglement between A and B is truly 
nonlocal, may be used for teleportation, and will violate 
Bell inequalities. 

On the other hand, the same example shows that the 
amount of entanglement depends on the definition of the 
modes. Namely, given a set of mode functions Fj and 
corresponding mode operators at we are free to define 
new mode operators and functions j' by 



j 

Fj' = Uj>jFj , 



(6) 



where Uj/j is an arbitrary unitary matrix, so that condi- 
tions (§) and (g) are still satisfied [Q. In particular, we 
may rewrite the state (g) as 



M = |1>A'|0>B' 

in terms of new modes A' , B' , where 
a\, = (4 + a^)/V2, 



(a 



(7) 



(8) 



There is no entanglement in the state (Q) between sys- 
tems A' and B'. For example, if A and B correspond 
to horizontal and vertical polarization (with all other 
quantum numbers being the same), then the state (g) 
describes a single diagonally polarized photon which one 
normally would not call entangled. 

It is important to note that, depending on the context, 
not all redefinitions of the form (|6|) may be relevant. En- 
tanglement between two particular modes and the very 
definitions of the modes are only useful if one can perform 
measurements on those modes. With this criterion, two 
spatially separated modes A and B cannot be redefined 
in a useful way if only local measurements are possible. If 
we have two sets of modes {A\ . . . A n } and {B\ . . . B m } 
located in two different spatial regions then within each 
set redefinitions are allowed. However, such "local" re- 
definitions, as expected, cannot change the nonlocal en- 
tanglement between the two sets of modes. In particular, 
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in most (if not all) experiments entanglement is created 
between different spatial modes. For example, in down- 
conversion experiments || the so-called signal and idler 
modes propagate in different directions and may be en- 
tangled. Similarly, a two-mode squeezed state |(| is pro- 
duced by splitting two squeezed light beams on a beam 
splitter M with the two output beams propagating in dif- 
ferent directions. In such cases there is a preferred set of 
modes, and for states with a definite number of photons 
one does get the correct amount of entanglement between 
those modes even if one treats photons as systems. 

Perhaps more importantly, one can view the redefini- 
tion as a unitary operation that one can apply actively 
in order to create (or destroy) entanglement. For exam- 
ple, a polarizing beam splitter may be used to turn the 
state (pp into the state rtq). In general, if all modes in the 
problem are spatialy different then all unitary operations 
of the form (||) can be performed with just linear optics 
(beam splitters and mirrors). Other cases not covered by 
polarization or spatial degrees of freedom are d iscus sed 
in the Examples Section (specifically, in Section II C). 

We will discuss several examples of pure states and give 
the maximum and minimum entanglement possible by 
varying over all possible unitary redefinitions/operations 
of the form (^J) . We denote these quantities by E max and 
-Emin- We only consider transformations that leave the 
number of modes the same (or at least that do not in- 
crease the number) . Otherwise, any state with more than 
a single photon can be transformed into one with an ar- 
bitrarily large amount of entanglement by introducing an 
arbitrarily large number of modes. We also restrict our- 
selves to bipartite entanglement. In the simple example of 
the state (§) it is easy to see that under these conditions 
E m in = and £ max = 1. 



II. EXAMPLES 

We first consider states where the total photon number 
is fixed and equal to 2. 



A. 1 20) + |02) 

Consider a state where we have two modes and two 
photons in a state 



|*) = (|2) A |0) fl + \0)a\2)b)/V2. 



(9) 



For example, the two modes might describe horizontal 
and vertical polarization (with the remaining quantum 
numbers the same for both modes). Writing this state, 
incorrectly, as 

l*> = (I 1)1 I) + I ~>l ~>)/V2 (wrong!), (10) 

we would conclude this state has one ebit of entanglement 
and we might even think that that is the only correct an- 
swer (Note here that this notation is usually meant to 



imply that the two kets refer to different spatial modes; 
in that case the notation is not wrong, just a little dan- 
gerous). However, on the one hand the state @ can be 
rewritten as a product state 



I*) 



|1)4<U) 



so that E mm = 0, by redefining 

a\, = (a\ + ia ] B 
a B' = i a A - ia B 



(11) 



(12) 



Thus, when the original modes A and B describe linear 
polarization (as above), the state (|^) is equivalent to a 
state containing one left-hand circularly polarized photon 
and one right-hand circularly polarized photon. On the 
other hand, we have = log 2 (3); we can also rewrite 

the state as 



|*> = (|1>A«|1>*" 



\0)a»\2) b „)/V3. (13) 



The unitary operation accomplishing this transformation 
is less obvious now, but one of them is 



{a\-. 



(14) 



where x = 1/3 + 2\/2«/3. One cannot get more than 
log 2 (3) ebits of entanglement between 2 single modes, 
because the reduced density matrix of any mode can be 
at most of rank 3. 



B. |0110) + 1 1001) 

Now consider two photons in 4 modes A. . . D. Suppose 
we have a state of the form 

|*) = (|0>|1)|1)|0) + |1)|0)|0)|1»/V2, (15) 

where we left out the bothersome subscripts ABCD to 
indicate the modes. By writing this state in dangerous 
notation as 



\*) = (\B)\C) + \A)\D))/V2, 



(16) 



we may conclude that the entanglement between the pho- 
tons is 1 ebit. Indeed, when we restrict ourselves to bipar- 
tite entanglement between two pairs of modes, it turns 
out that Emin = 1- But in this case E maK = 2. Namely, 
we can rewrite the state jl5| ) as 



|*) = (|0)|1)|1)|0) + |1)|0)|0)|1) 
+ |0)|0)[|0)|2) + |2)|0)]/V2 
- [|0)|2) + |2)|0)]|0)|0)/V2)/2, 



(17) 



which manifestly possesses 2 ebits of entanglement be- 
tween the first pair of modes A 1 , B' and the second pair 
C',D'. The corresponding transformation is 
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x + a\ — X-O* 



Di 



X-a' A + x + a* D , 



(18) 



with x± = y 1/2 ± \/2/4. Since p^'S' can be at most of 
rank 4, 2 ebits of entanglement is the maximum possible. 

The state (|l5]) possesses one ebit of bipartite entangle- 
ment not only between any two pairs of modes (chosen 
from ABCD) but between any one single mode and the 
remaining three. As a curiosity we note that under uni- 
tary redefinitions the minimum bipartite entanglement 
between one mode and three others is no longer 1 ebit 
but less, E min = 2 - 31og 2 (3)/4 w 0.8113. In this mini- 
mum, the reduced density matrix of the single mode is 



p=i[3|0)(0| + |2)<2|]. 



The maximum possible entanglement is E, 
which was found numerically. 



(19) 

,„.* « 1-3002, 



C. 100001) + 1010010} + 1001100} 



in ||, where the entanglement was called "double en- 
tanglement" , referring to the fact that the entanglement 
appears to be equal to 2 ebits. 

Let us conclude this subsection by discussing briefly 
physical implementations of a state like (pp|). The six 
modes A. . . F could, of course, correspond to different 
spatial modes, or to two orthogonally polarized sets of 
three different spatial modes. But a recent experiment 
H suggests one may alternatively make use of a different 
degree of freedom, namely the transverse mode profile: 
Within the paraxial approximation, a beam propagating 
in a given (z) direction, with fixed polarization and fixed 
frequency, still has an infinity of degrees of freedom left. 
For instance, one complete set of paraxial mode functions 
is given by the Hermite-Gaussian (HG) modes, another 
by the Laguerre-Gaussian (LG) modes p0| . The latter 
modes may be characterized by a quantum number L z , 
describing orbital angular momentum of light 1 1 1 . It is 
not straightforward to implement arbitrary unitary op- 
erations of the form @ in this case. However, specific 
transformations are known: the conversion of HG modes 
into LG modes of the same order and vice versa with 
astigmatic lenses was discussed in Refs. Uj, and a larger 
set of unitary operations is considered in Ref. |l3[ ] . En- 
tangled states using the transverse degrees of freedom 
were discussed in [pTf . 



Continuing the series with two photons in 6 modes 
where the state has the symmetric form 

|*) = (|0)|0)|1)|1)|0)|0) + |0)|1)|0)|0>|1)|0) 

+ |1)|0)|0)|0)|0)|1))/V3, (20) 

generalizing (|l5|), we might conclude from writing the 
state as 

\V) = (\A)\F) + \B)\E) + \C)\D))/V3, (21) 

that we have log 2 (3) ebits of entanglement, since each 
photon "has a three-dimensional Hilbert space attached 
to it" , to use the incorrect picture once again. 

However, restricting ourselves to bipartite entangle- 
ment between two triplets of modes, we have E min = 1 
and E max — log 2 (5). This is the direct generalization 
of the result of the preceding subsection. Here, log 2 (5) 
ebits is the maximum possible because the reduced den- 
sity matrix pA',B',c can be at most of rank 5. 

We in fact conjecture that for states with 2 photons 
in 2N modes of the form (|2l|), or (po|), (which appear 
to have log 2 (iV) ebits of entanglement) we always have 
-^min = 1 an< 4 -^max = log 2 (JV+2), if we consider only en- 
tanglement between two sets of N modes. Note that the 
reduced density matrix of N modes in a state with ex- 
actly N photons can be at most of rank N+2: there are N 
different ways in which 1 photon can be distributed over 
N modes, and there are 2 ways either set of N modes 
can have no photons. The case N = 8 was considered 



D. 0220) + |2002) - |1111) 

A 4-photon state of the form 
|*) = (|0)|2)|2)|0) + |2)|0)|0)|2) - |1)|1)|1>|1»/V3, (22) 

features in a recent experiment [ p"5| . There are two pre- 
ferred pairs of modes in the experiment (the first two 
and the last two) and so it is correct to infer that the 
state ( |22] ) possesses log 2 (3) ebits of entanglement be- 
tween those two pairs of modes. Restricting ourselves 
to bipartite entanglement between pairs of modes, this 
is in fact equal to the minimum, E m - ln = log 2 (3). By 
numerically searching for the maximum we found that 
-E max ~ 2.9798, which is obtained when the reduced den- 
sity matrix of a pair of modes is of the largest possible 
rank 9 |ll|, but not maximally mixed (log 2 (9) w 3.17). 

Since in the experiment of Ref. fl5[ ] the 4 modes in- 
volved are orthogonally polarized signal and idler modes, 
only ordinary and polarizing beam splitters are needed 
to transform the original state with log 2 (3) « 1.585 ebits 
into one with almost 3 (in fact, E max ) ebits of entangle- 
ment. 



E. 1 00) + |11) 

Our final example is a superposition of states with dif- 
ferent photon numbers, 
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|*> = (|0>|0> + |1>|1»/V2. 



(23) 
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Again, this example shows that the entanglement is not 
between photons since the presence of the vacuum is 
here an essential part of the entanglement. Written in 
the form (p3| ) we would conclude we have 1 ebit of en- 
tanglement. However, for the minimum amount of en- 
tanglement we find E min w 0.3546. More precisely, we 
get E min = - J2 k X k log 2 A fc with Ai, 2 = 1/2 ± ^3/4 the 
eigenvalues of the reduced density matrix of a single mode 
if we rewrite the state as 

l*> = (|0)|0) + (|0>|2) + |2>|0))/^)/V2, (24) 

which is obtained by the same transformation (||) as in 
Subsection II A . The maximum amount of entanglement 
is, surprisingly perhaps, just a tiny little bit over 1 ebit, 
namely E max w 1.0071, as was found numerically. 



III. CONCLUSIONS 

There are two differences between a second-quantized 
theory, such as Quantum Electrodynamics, and a first- 
quantized theory, such as standard non-relativistic quan- 
tum mechanics of material particles. One is that parti- 
cles, photons in particular, are not systems in a second- 
quantized theory. The other difference relevant to 
quantum-information processing is that the systems, the 
EM field modes in particular, can be redefined arbitrarily, 
at least in principle. But one does not have any freedom 
in (re)defining ions in an ion trap or atoms inside a cav- 
ity. (Of course, in a second-quantized theory of matter 
particles this would be different, but we live in a world 
where non-relativistic quantum theory is an extremely 
good and convenient approximation for atoms and ions.) 

It turns out that the amount of entanglement present 
in a given state depends on how one defines one's sys- 
tems. The unitary transformations (^) may be viewed 
as, indeed, merely redefinitions of field modes. In prac- 
tice, however, there is a more useful point of view. The 
transformation (0) may be seen as a unitary operation on 
modes that can actually be performed by linear optics el- 
ements such as beam splitters, polarizing beam splitters, 
or in more complicated cases by astigmatic lenses. This 
way the entanglement in a given state may often be in- 
creased beyond the amount that is "obviously" present 
in the state, as was shown in several examples. 

Finally, let us point out some relations with previ- 
ous work. The conclusions reached in the present paper 
agree, in some indirect sense, with those in two quite dif- 
ferent papers. First, in Ref. |l7j it is shown that certain 
aspects of local entanglement, that is, entanglement be- 
tween different modes that are all in the same location, 
can be simulated with classical light beams. Second, in 
Ref. Jl8| it is shown that the amount of entanglement 
between the spatial and spin degrees of freedom of a rel- 
ativistic spin-1/2 particle depends on the reference frame 
of the observer. 
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